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Abstract Observations carried out with the Magneto- Optical Filter at Two 
Heights (MOTH) experiment show upward-traveling wave packets in magnetic 
regions with frequencies below the acoustic cut-off. We demonstrate that the 
frequency dependence of the observed travel times, i.e. the dispersion relation 
shows significant differences in magnetic and non-magnetic regions. More impor- 
tantly, at and above the layer where the Alfven speed equals the sound speed 
we do not see the dispersion relation of the slow acoustic mode with a lowered 
cut-off frequency. Our comparisons with theoretical dispersion relations do not 
suggest this is the slow acoustic wave type for the upward low-frequency wave. 
From this we speculate that partial mode conversion from the fast acoustic to 
the fast magnetic wave might take place. 

Keywords: hclioscismology; wave propagation; magnetoacoustic waves; radia- 
tive transfer 

1. Introduction 

Recently, the observation and analysis of upward traveling waves in the solar 
atmosphere have become of great interest, as they are thought to provide a 
considerable amount to the heating of the chromosphere. Commonly it was 
understood that only waves with frequencies above the acoustic cut-off fre- 
quency (luq) could travel freely in the solar atmosphere, whereas waves with 
lower frequencies are trapped inside the acoustic cavity of the Sun. However, 
recent studies by several researchers have detected upward traveling waves with 
frequencies lu < luq. 

First, from observations with the Transition Region and Coronal Explorer 
(TRACE; Handy et al, f999) Dc Pontieu, Erdelyi, and James (2004) and Mcin- 
tosh and Jefferies (2006) show that the presence of an inclined solar mag- 
netic field allows low-frequency waves to propagate into the solar chromosphere. 
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The same effect was shown by Jefferies et al. (2006) with data taken with the 
Magneto- Optical filter at Two Heights (MOTH; Finsterle et al. 2004a) experi- 
ment. Furthermore, Vecchio et al. (2007) found low-frequency upward-traveling 
waves in the velocity signal measured simultaneously in Fei 7090 A and Can 
8542 A with the Interferometric BIdimensional Spectrometer (IBIS; Cavallini 
2006). 

What is the nature of the observed low- frequency waves, i.e. are these Alfven 
waves, slow acoustic, or fast magnetic waves? De Pontieu, Erdelyi, and James 
(2004), Mcintosh and Jefferies (2006), and Jefferies et al. (2006) consider the low- 
frequency waves in areas with inclined magnetic field to be the ficld-guidcd slow 
acoustic waves mode (see Section[2] for details) leaking through "portals" where 
the inclined magnetic field lowers the acoustic cut-off. Numerical simulations, 
carried out by Steiner et al. (2007) and Rosenthal et al. (2002) indicate the pos- 
sibility of mode conversion and transmission, which goes in line with theoretical 
calculations by Cally (2005, 2006, 2007). From MOTH observations, Finsterle 
et al. (2004b) show that the "magnetic canopy" (i.e. the (3 w 1 layer, with 
/? = 87rp g /B 2 , and p g being the gas pressure and B the magnetic flux density) 
has a strong influence on the propagation behavior of magnctoacoustic-gravity 
waves. 

2. Wave Modes 

The different wave modes in the solar atmosphere are very well understood from 
the theoretical point of view, however the determination of the wave modes from 
the observation remains still a difficult task. The actual wave type depends on the 
local sound speed (c), Alfven speed (a), pressure (p), temperature (T), magnetic 
flux density (B) and field inclination (9). Generally, the following types of waves 
can be distinguished (for details see e.g. Rosenthal et al., 2002; Bogdan et al., 
2002, 2003; Khomenko and Collados, 2006; Cally, 2005, 2006, 2007; Schunker 
and Cally, 2006). In high plasma-/? regions (a < c), the fast mode is a longitu- 
dinal acoustic wave. In low plasma-/? (a > c) the slow mode is a longitudinal 
acoustic wave propagating along the magnetic-field lines. Here, the fast mode 
is the magnetic mode, which, in the case of relatively weak magnetic fields, is 
nearly aligned with the wave vector irrespective of the field direction. However 
for high magnetic-field strengths, the fast-mode displacement is perpendicular 
to B irrespective of the propagation direction (Rosenthal et al., 2002). Finally, 
the Alfven wave is always transverse to both the magnetic field and the wave 
vector. 

Stein (1971), Rosenthal et al. (2002), Bogdan et al. (2002, 2003), and Steiner 
et al. (2007) showed from 2D numerical magneto-hydrodynamics (MHD) sim- 
ulation that the fast and slow waves undergo coupling at the layer where the 
plasma-/? is of the order of unity. Similarly, Cally (2005, 2007) suggested that 
mode transmission and mode conversion take place at the layer where the sound 
speed equals the Alfven speed. In Cally's terms, mode transmission is the change 
of a fast acoustic (magnetic) wave to a slow acoustic (magnetic) wave. Mode 
conversion is the change from a fast (slow) acoustic to a fast (slow) magnetic 
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Figure 1. Phase travel time map (455 Mm X 455 Mm at disk center) for waves with frequency 
o;=3.0mHz (a) and the concurrent MDI magnetogram (b). The red and yellow pixels indicate 
upward traveling waves with phase travel times up to 22 seconds, coinciding with magnetic 
regions in the MDI magnetogram, whereas t p h Rj for the non-magnetic areas. 



wave. In the following we use Cally's notation of mode conversion and trans- 
mission. It is important to note that (3 = (2/^)c 2 /a 2 « 1.2c 2 /a 2 . This means 
that the /3 = 1-laycr and the layer where a = c are contiguous. We use the term 
conversion layer for the layer where a = c. This condition reflects the matching 
of the phase velocities of the slow and fast waves, which is the fundamental 
physical precondition for the waves to couple (Cally, 2005). In order to identify 
the wave type from the observation "it is imperative to determine whether the 
spectral line or continuum is formed in low or high-/? plasma" (Bogdan et ai, 
2003). This is pursued in the present paper with respect to the conversion layer. 



3. Data 

The MOTH experiment observed simultaneous Doppler shifts in the K I 7699 A 
and Nai 5890 A lines during a South Pole campaign in the austral Summer 
2002-03 (Finstcrlc et al, 2004a) ofthcfull solar disk with a resolution of3.7Mm 
pixel -1 . Here we make use of the (455 Mm x 455 Mm 17.8 hour) data cube taken 
at disk center starting on 2003 January 20, 00:59 UT. Figure 1, panel (a) shows 
the phase travel time £ p h between the K I 7699 A and Na I 5890 A observing layers 
for waves with frequency cj=3.0mHz. In non-magnetic regions the mean phase 
travel time is i p h ~ (green-blue areas), while the red and yellow pixels indicate 
upward traveling waves in regions that coincide with enhanced magnetic flux 
density in the concurrent MDI magnetogram (panel b). Under the assumption 
of a significant amount of inclined field lines in the magnetic areas, this agrees 
well with the findings by De Pontieu, Erdelyi and James (2004), as well as 
Mcintosh and Jefferies (2006). 



Haberreiterl2_arxiv.tex; 28 January 2010; 1:21; p. 3 



1 



Haberreiter and Finsterle 



Table 1. Formation height (km) with respect to T5000 = 1 calculated for 
the quiet Sun and plage for disk center. Also listed are the mean sound 
speed (c) at the formation heights and the cut-off frequency (uiq) derived 
from the atmosphere structures. 



area ZNa (km) 2k (km) Az (km) c (kms 1 ) ujq (mHz) 



Quiet Sun 530 160 370 6.5 5.6 

Plage 490 270 220 7.2 5.0 



4. Conversion Layer and Formation Heights 

Purely acoustic waves can not propagate below the acoustic cut-off frequency. 
Following the concept of mode transmission and conversion (Cally, 2007) im- 
plies that low-frequency upward-traveling waves only exist in regions where the 
conversion layer is below the formation heights of the spectral lines employed in 
the MOTH experiment. Otherwise the low-frequency wave would be evanescent. 
The reason is that only if the conversion layer is below the height of observation 
can the upward-propagating wave have coupled to the magnetic field. This goes 
in line with the concept of portals by Jefferies et al. (2006). To further test this 
hypothesis we use magnetic-field extrapolations (Mcintosh et al., 2001) of the 
MDI magnetogram in Figure lb to estimate the plasma-/? from the plage model 
atmosphere structure by Fontenla et al. (2006). The field extrapolations are car- 
ried out with respect to the mean observational height of Ni 1 6768 A in the four 
MDI-filters employed to determine the magnetic-flux density. Furthermore, from 
radiative transfer calculations carried out with COSI (Haberreiter, Schmutz, and 
Kosovichev, 2008; Haberreiter, Schmutz, and Hubeny, 2008) we determine the 
formation height of K 1 7699 A and Na 1 5890 A at disk center for the quiet Sun 
and plage. Note that the formation heights change for different view angles or 
positions on the solar disk. Here we focus only on plage and leave the sunspots 
for a later study. 



From t\ = 1 and the position of the Magneto-Optical Filters (MOFs) in the 
MOTH instrument we derive the mean observing height of K 1 7699 A and Na 1 
5890 A, given in Table[T] For a detailed discussion of the MOTH observing heights 
see Haberreiter, Finstele, and Jefferies (2007). Note that we do not account for 
any magnetic effects on the formation of the line. This will be pursued in future 
work. Figure[2] shows the contours of the calculated conversion layer below the 
Nal 5890 A observing height for plage. The red contours indicate where the 
conversion layer coincides with the observational height of Nal 5890 A (490km). 
Toward the center of the active regions the conversion layer drops and reaches 
heights lower than the observational height of Ki 7699 A, zk = 270 km (green 
to blue contours). The striking congruence of the areas where waves propagate 
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Figure 2. Heights of the conversion layer (a = c) where it is below the Nal 5890 A formation 
height in plage with respect to /i(t5ooo = 1)- The areas where the conversion layer is below 
the Na I 5890 A formation height coincide with areas of higher phase travel time in Figurc[l] 
White indicates that the conversion layer is above the formation height of Na. 



(FigurcQJi) and where the conversion layer cuts through the MOTH observing 
layers, in particular the areas where z a=c < zk supports the assumption by Jef- 
fcrics et al. (2006) that the low-frequency waves propagate only in low-/3 regions. 
Moreover, this consistency suggests that the formation heights calculated with 
COSI from the ID atmosphere structures are correct. 



5. Dispersion Relations 

Low-frequency traveling waves in the solar atmosphere are generally explained 
by a lower cut-off frequency due to an inclined magnetic field. According to 
Campos (1987) the reason is that at high altitudes (i.e. low-/?) the acoustic- 
gravity waves must follow the magnetic-field lines, increasing the scale height by 
the factor (cos0) -1 , with 8 being the field inclination to the vertical. Formally, 
this can also be described with a change of the effective gravitational acceleration 
<?eff = go cos 8 (Dc Pontieu, Erdelyi, and James, 2004). A lowering of the acoustic 
cut-off frequency, however, implies that the propagating wave should still show 
the acoustic dispersion relation. To test this hypothesis we compare the group 
and phase travel times derived from the MOTH data with theoretical values. 
The fitting of phase travel time (£ p h) and the group travel time (i gr ) is described 
in detail by Finsterle et al. (2004b) and briefly outlined here. In order to study 
the frequency dependence of the wave propagation, the K and Na data cube is 
first frequency-filtered in the time domain with a Gaussian of width 0.57 mHz. 
Note that, as wc are interested in the study of the dispersion of propagating 
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wave-packets, the finite width of the frequency filter is necessary to obtain the 
superposition of waves with different frequencies and at the same time be able 
to localize the wave packet. The width of 0.57 mHz was chosen, as it allows 
the localization of a wave-packet with a reasonably well defined frequency. For 
our analysis we apply a frequency filter between 2.5 and 7.5 mHz in steps of 
0.285 mHz, indicated by the crosses in Figure^. The frequency-filtered time 
series are then cross correlated. Finally, the resulting cross-correlation functions 
are fitted by least-square fits to 



, , A 2 8uj 
y(t) = ~^=- exp 

VSTT 



5u 2 (t-t gr ) 2 



cos[w(t - i p h)] (1) 



where i pn , i gr , the amplitude (A), the central frequency (w), and the width of 
the Gaussian filter (Suj) are the fitted variables. In this analysis the filter width 
dui is allowed to vary as the effective width of the filter is also influenced by the 
shape of the underlying spectrum. Nevertheless, the fitted widths are typically 
very close to the given values. The fitted Garbor wavelet function allows us 
to directly determine the group and phase travel time. We note that a cross- 
spectrum analysis would be an equivalent procedure to derive these parameters. 

In Figure^ the fits for t gr (thin lines) and i p h (thick lines) are shown for 
different height intervals of the a = c-layer (z a — c ) with respect to the formation 
height of Na (zNa) and K (zk)'- 

i) z a=c > ZNa (i-e. a < c, blue lines), 

ii) ZNa > z a=c > zk (i.e. a sa c, red lines), and 
Hi) zk > z a=c (i.e. a > c, black line). 

We find the dispersion relation of acoustic-gravity waves (u>o ~ 5.4 mHz) for the 
heights where a < c (blue lines), showing that the criterion a < c is suitable 
to identify quiet- Sun conditions and that it allows to reproduce the results by 
Finsterle et al. (2004b). For the regime where a w c we find traveling waves 
at low frequencies (w < 5 mHz) with a mean phase travel time of 5.2 seconds 
(thick red line), and for a > c a mean phase travel time of 9.1 seconds (black 
line), i.e. with increasing magnetic- flux density the phase travel time increases. 
For this analysis we consider areas with magnetic-flux density up to B =500 G. 
The mean magnetic-flux density for the different height intervals are B=8 G 
(z a=c > z Na ), 5=109 G (z Na > z a=c > z K ), and 5=250 G (z a=c < z K ). Note 
that as the filling factor for plage regions is smaller than unity, at some areas 
below the instrumental resolution the "real" magnetic flux density is expected to 
have higher values than B. The corresponding la-error bars of the travel time 
fits are shown in Figurc[3]3,c. Moreover, with our approach cannot determine 
whether the signal observed in one pixel is a superposition of different waves. 
Due to the fact that the filling factor can be less than one this cannot be ruled 
out. However, from our analysis it is clear that for the a c regime the data do 
to reveal the predominant signal of slow magneto-acoustic gravity waves with a 
lowered cut-off frequency. 

The group travel time for the a > c-regime is not shown due to the increased 
noise in the fitting of the group travel time combined with poor statistics due to 
the sparce number of pixels with higher magnetic-flux density. In the following 
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we discribe the comparison of these results with theoretical dispersion relations. 
5.1. Acoustic Wave (a <C c) 

For the fast acoustic (fa) wave we use the dispersion relation 

LO 2 = C 2 k 2 +L0l, (2) 

going back to Bel and Leroy (1977), with c being the sound speed, (k = a + if3) 
the complex wave number and ujq = c/4H is the cut-off frequency for an idealized 
isothermal atmosphere. Generally, the cut-off frequency for a ID non-isothermal 
atmosphere is u> 2 = c 2 / \AH 2 )(1 — 2dH/dz). As dH/dz is small in the layers 
where Nai 5890 A and Ki 7699 A are formed, we consider the approach by Bel 
and Leroy (1977) justified. Note, that Equation @ can also be derived from the 
dispersion relation given by Schunker and Cally (2006), Equation (12), for the 
limit of the non- magnetic case, i.e. a — > 0. 

We are aware that the assumption of vertically upward traveling waves prop- 
agating long inclined field lines is an approximation that is of course not always 
the case on the Sun. However, for the MOTH instrument the observation of 
the horizontal component of the wave is limited to small inclination angles. 
Let us discuss this briefly with the following example. Consider a wave with 
v =3.5 mHz. With a sound speed of 7kms _1 we get a wavelength A =2000 km. 
The horizontal sampling rate (S) of the MOTH data is Rj3700km pixel -1 . The 
maximum angle (6> max ) (6 being the angle of the propagation vector and the 
vertical) under which the horizontal component of the wave can still be resolved 
with MOTH is: sin((9 max ) = A/(2S) « 0.27. From this we get max ps 16°. Note 
that for higher frequencies 8 max is smaller, i.e. for 7 mHz it is rj 8°. As these 
angles are small, we consider the ID approach appropriate. 

From Equation @ we get k = ±^/lu 2 — loq/c and the phase velocity: 



Wph.fa — 




Furthermore, using u = ±y c?k 2 + ujq we get the group velocity: 



^gr,fa — 




For u! > wo, k = a and /3 = 0, while for u < u>o, k is imaginary (fc = i(3), meaning 
it denotes the evanescent wave, i.e. the wave varies sinusoidally in time, but 
exponentially in space. In the latter case, for a single component wave, a = 0, 
u p h = oo, and v gI = 0. With i p h = Az/v p h and t gr = Az/v gI it follows that for 
the evanescent wave we expect i p h = and t gr = oo. However, as it is impossible 
to observe an infinite group travel time, we set the theoretical group travel times 
for the evanescent wave to zero. This is justified by the fact that in the case of 
the evanescent wave the whole solar atmosphere oscillates in phase, i.e. there 
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Figure 3. Panel a: Group and phase travel time (thin and thick lines) derived from the MOTH 
observations as a function of frequency for different height intervals: a < c (blue lines), a fa c 
(red lines), and a > c (black line). Panel b: For a < c-regime we find the dispersion relation as 
expected for an acoustic-gravity wave with a cut-off frequency. Also shown are the theoretical 
travel times (t gr and t p h) for an acoustic gravity wave with a cut-off frequency luq = 5.6 mHz, 
which agrees well with the observed travel times. Panel c: For the regime where a ss c we find 
positive travel times below the acoustic cut-off. Also shown are theoretical dispersion relation 
based on a lowered cut-off frequency uj e g = wrjcosf, where 8 = 30° and uio = 5.0 mHz. 
The observed travel times are significantly different from the theoretical values, indicating a 
different wave type. Panel b and c also show the l-cr error bars for the fitted travel times. Note 
the different scale of the y-axes. 
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is no detectable phase lag between different layers in the atmosphere. As we 
will see later, this approach agrees well with the observed dispersion relation 
in non-magnetic regions. Note that under the assumption of an infinitely wide 
region where a wave is evanescent, the wave does not transmit energy vertically 
(e.g. Lcibacher and Stein, 1981). However, if this regime is finite, tunneling can 
occur, in which case the wave does carry energy and there is phase propagation 
vertically. 

The acoustic cut-off frequency uq = jg/Aitc is calculated from the latest 
solar atmosphere structures for the quiet Sun (Fontenla, Balasubramaniam, and 
Harder, 2007) and plage (Fontenla et ai, 2006). Furthermore, 7=5/3 is the ratio 
of the specific heats, £/=274ms -2 the gravitational acceleration, c = V7p/p the 
sound speed, p and p the depth dependent pressure and density. Note that this 
is only one of several possible representations of the acoustic cut-off frequency 
in the quiet Sun (Schmitz and Fleck, 1998). 

5.2. Magnetic Waves (a ^> c) 

Let us now discuss the dispersion relation of upward traveling wave modes that 
can occur after mode transmission or conversion, which are the field guided slow 
acoustic mode and the fast magnetic mode. In the present study we focus on 
the question of whether we can identify the slow acoustic mode. This mode 
is understood to be "leaking" through magnetic regions by the lowering of the 
acoustic cut-off frequency in magnetic regions (see Bel and Leroy, 1977; Campos, 
1987; De Pontieu, Erdelyi and James, 2004; Cally, 2006). Here we account for 
the lowering of the acoustic cut-off frequency following Campos (1987, Equation 
209a) and Dc Pontieu, Erdelyi and James (2004): 

u 2 = c 2 k 2 +ul Sl (5) 

with cj c ff = ujq cos 8, and 9 being the field inclination derived from the magnetic- 
field extrapolation. Note that this dispersion relation can again be derived from 
the dispersion relation given by Schunker and Cally (2006), Equation (12) for 
the limit of the magnetic case, i.e. a 3> c. Following the above outline for the 
acoustic wave, we get the phase and group velocity v p h, S a and v griS a for the field 
guided slow acoustic (sa) wave mode as: 



^ph,sa — 




tWa = C^l-(^) . (7) 

6. Results 

In the following we compare the travel times derived from the observations 
with the theoretical travel times as outlined above. For consistency with the 
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observations, we apply a running frequency filter to the theoretical travel times. 
This filter is the same Gaussian filter as used for the frequency filtering of the 
MOTH data (6w = 0.57mHz). FigurejSbjC compares the observed travel times 
with the group and phase travel times, calculated as i p h(w) = Az/v p h(oj), and 
ig r (cj) = Az/Vgr(oj) for different regimes in the solar atmosphere. Panel (b) 
compares the observed travel times for a < c with the theoretical values calcu- 
lated with A^=370km, c =6.5 kms— 1, ojq = 5.6 mHz (see Table 1) and no field 
inclination (8 = 0°). The theoretical curve represents the observed dispersion 
relation quite well. For the intermediate-/? regime (a « c, panel c) we calculate 
the dispersion relation using Az=220km, c —7.2 km, and ujq = 5.0 mHz for 
plage, and a lowered cut-off frequency with 9 — 30° (dashed lines) . We find that 
the fitted t gI shows a substantially different behavior than the theoretical t gI . 
In particular, a cut-off frequency is not visible in the travel time fits from the 
observation. Therefore, we conclude that, while passing through the magnetic 
canopy, the acoustic nature of the wave is no longer detectable. 

7. Discussion 

Let us discuss the significant change of the dispersion relation seen in Figurc|31 
panels (b) and (c). As the observation heights of Na and K are very close to the 
conversion layer, one might argue that we do not detect a pure slow acoustic 
wave with a lowered cut-off frequency, but a mixture of wave types, i.e. the slow 
acoustic and the fast magnetic. However, the clear correlation signal between 
both observing layers does not support this point. Therefore, it is striking that 
the acoustic dispersion relations seen in the high-/? regime vanishes immediately 
when the observing heights are close to the conversion layer. It is important to 
note here that the cross correlation yields robust results, so the wave packets 
lead to a clearly correlated signal in both observing layers. 

Furthermore, as the gas and magnetic pressure are of approximately equal 
amplitude, i.e. the magnetic field is weak, it is possible that the acoustic wave 
coming from below not only causes longitudinal pressure and velocity fluctua- 
tions along the field lines, but also distorts the magnetic-field lines along the 
wave vector (see Rosenthal et at, 2002). The combination of these processes 
might lead to a loss of the acoustic signature in the observations. 

One might expect to see the pure slow acoustic wave if we could discriminate 
waves propagating upward under an angle parallel to the magnetic field. In this 
case the acoustic wave can not distort the magnetic-field lines. Holography, as 
carried out by Schunker et al. (2008) for vertically upward traveling waves would 
be the right tool to select the appropriate incoming waves. Using this technique 
we might be able to analyze the dispersion relation of an 'ideal' slow acoustic 
wave. 

The seemingly slight decrease of the travel times for the magnetic case (Fig- 
ure[3] dashed lines) for frequencies higher than 5 mHz also needs further investi- 
gation. Muglach, Hofmann, and Staude (2005) conclude from MDI and TRACE 
observations that if the conversion layer is below the formation height of the 
TRACE intensity, the intensity power in the layers above decreases, which they 
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associate with wave reflection. Numerical calculations carried out by Rosenthal 
et al. (2002) show that wave reflection of high frequency waves (in their case 
uj = 42 mHz) occurs at the (3 = 1 layer. In the regions above the reflecting 
surface the wave is still present as an evanescent tail. This could explain the 
vanishing travel times in the high frequency regime. However, only detailed 3D 
MHD simulations including realistic density and temperature stratifications can 
help to fully understand these phenomena. 

8. Conclusion 

From MOTH observation we show that in magnetic regions low-frequency waves 
can escape the solar acoustic cavity and propagate upwards between the K I 
7699 A and Nai 5890 A observing layers. Furthermore, evidence is given that 
the observed low-frequency waves originate in deep layers as acoustic (c 3> a) 
and clearly change their characteristics while passing through the conversion 
layer into the low- (8 regime. In particular, we do not see any evidence for a 
lowered cut-off frequency in low-/? regions. The reason for this is currently not 
explained and needs further investigation. 

It is however clear that in order to distinguish various wave types in the solar 
atmosphere, for a detailed analysis the observational heights have to be sepa- 
rated into the very high and low-/3-regimes. This is essential for new instrument 
designs. Only by simultaneously cutting through very high and low-/? regimes 
can we analyze the characteristics of the waves in detail. As the different wave 
types dissipate their energy in different heights of the solar atmosphere, it is 
essential to determine to what extent the waves convert and transmit into the 
slow magnetic and fast acoustic wave. 
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Appendix 

In the following, we show how the dispersion relations given in Equations. Q 
and ([5]) can be derived from Equation (12) in Schunker and Cally (2006). Let us 
first consider the acoustic wave, i.e. a — > 0. Under the assumption of vertically 
upward traveling waves, i.e. k x —0, their equation simplifies to: 

uj 4 - (a 2 + <?)k 2 u? + a 2 c 2 k 4 cos 2 6 - (lo 2 - a 2 k 2 cos 2 6>)w 2 = 0. (8) 

Using the asymptotic case of a = 0, Equation ([8]) further simplifies to 



uj 4 - c 2 k 2 ui 2 - uj 2 ^ 2 = 0. (9) 

After dividing with w 2 , this then leads easily to Equation @ . Now let us consider 
the asymptotic magnetic case with a >> c or 1 /a 2 — > 0. Multiplying Equation (J5]) 
with 1/a 2 we get 

- k 2 uj 2 + c 2 k 4 cos 2 9 + k 2 cos 2 6uj 2 = 0. (10) 
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Dividing by k 2 and using cos 2 9cu 2 = lu 2 s , and assuming vertically propagating 
waves leads to Equation ([5]). 
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